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Abstract.  Transient  motions  of  a  vapor-gas  mixture  due  to  the  evaporation  and  condensation  processes  from  or  onto  the  plane 
condensed  phase  with  a  temperature  field  as  its  internal  structure  have  been  studied  numerically  based  on  the  Boltzmann 
equation  of  BGK  type  for  a  binary  gas  mixture.  The  condition  of  the  continuity  of  the  energy  flow  across  the  interface  surface 
has  to  be  imposed  in  addition  to  the  conditions  at  the  surface  of  the  condensed  phase  without  the  internal  structure  because  the 
temperature  of  the  surface  in  this  case  becomes  an  unknown  parameter  to  be  determined  as  part  of  the  solution.  The  transient 
flow  fields  of  the  mixture  occur  owing  to  the  phase  change  processes  caused  by  the  continuous  change  in  temperature  of  the 
surface  of  the  condensed  phase.  The  effects  of  the  internal  structure  of  the  condensed  phase  coupled  with  the  presence  of 
a  noncondensable  gas,  which  is  small  in  amount  here,  on  the  transition  process  of  the  transient  flow  fields  and,  hence,  the 
formation  and  the  propagation  of  the  shock  waves  and  the  contact  regions  may  clearly  be  recognized.  Some  of  the  features  of 
the  present  flows  are  that  1)  the  maxima  of  the  mass  and  energy  flows  may  exist  at  a  certain  value  of  the  latent  heat  parameter, 
the  phenomenon  of  which  may  be  ascribed  to  the  coupling  effects  of  the  latent  heat  parameter  and  the  existence  of  the  internal 
structure  of  the  condensed  phase;  2)  the  negative  mass  flow  appears  in  a  short  period  of  time  during  the  transitional  state  of 
the  flow  fields  when  the  latent  heat  parameter  is  extremely  small.  However,  it  disappears  soon  as  time  goes  on.  These  two 
phenomena,  which  are  affected  here  by  the  presence  of  a  noncondensable  gas,  are  also  found  in  a  pure- vapor  case. 


INTRODUCTION 

Phase-change  phenomena  and  flow  problems  associated  with  these  processes  are  not  only  of  theoretical  interest  in  gas 
dynamics  but  also  of  practical  importance  because  they  do  have  wide  applications  in  various  thermo-fluid  engineering 
fields.  Such  problems,  however,  are  the  ones  to  which  the  ordinary  continuum-based  fluid  dynamics  is  not  directly 
applicable  because  of  the  existence  of  a  nonequilibrium  region  called  the  Knudsen  layer  in  the  close  vicinity  of 
the  condensed  phase.  The  existence  of  this  region  can  never  be  neglected  in  any  problems  of  such  kind  even  in 
the  continuum  limit.  The  analyses  of  these  problems,  therefore,  must  necessarily  be  based  on  the  kinetic  equation 
because  it  is  this  nonequilibrium  region  that  is  responsible  for  the  phase-change  processes  to  occur.  However,  the 
kinetic  equation  is  complex  and  difficult  to  deal  with  and  the  analyses  based  on  this  equation,  regardless  of  whether 
they  are  theoretical  or  numerical,  impose  us  a  great  burden.  Moreover,  when  the  gas  phase  is  a  mixture  consisting 
of  a  condensable  and  noncondensable  gases,  the  problems  become  more  complex.  In  spite  of  this,  a  large  number 
of  evaporation  and  condensation  problems  have  been  studied  so  far  not  only  for  pure- vapor  cases  but  also  for  gas 
mixture  cases  (see  e.g.,  the  proceedings  of  the  Rarefied  Gas  Dynamics  Symposia)  based  on  the  kinetic  equation  and 
a  fairly  large  amount  of  the  useful  information  and  knowledge  has  been  obtained.  Among  others,  it  would  be  worth 
noting  from  the  practical  point  of  view  that  the  noncondensable  gas,  especially  when  it  is  of  small  amount,  is  driven 
toward  the  interface  and  accumulate  near  there,  forming  an  extremely  large  spacial  gradient  of  the  density  leading 
to  the  great  hindrance  of  the  phase  change  proesses  [1].  However,  the  analyses  so  far  are  all  confined  to  systems  in 
which  the  condensed  phases  have  no  internal  structures  except  a  recent  work  by  Onishi  and  Yamada  [2],  who  treated 
a  half-space  problem  of  a  pure- vapor  case  based  on  the  Boltzmann  equation  of  BGK  type  [3].  They  considered  a 
temperature  field  as  the  internal  structure  of  the  condensed  phase  and  studied  its  effects  on  the  flow  fields  of  the  vapor. 
The  consideration  of  the  internal  structure  within  the  condensed  phase  may  give  us  practically  important  information 
on  flow  fields  within  heat  and  energy  exchange  or  transfer  systems  such  as  heat  pipes.  Especially  for  mixture  cases 
containing  a  small  amount  of  noncondensable  gas,  the  effects  of  the  internal  structure  coupled  with  the  presence  of  the 
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noncondensable  gas  would  become  very  important  and  should  be  studied  because  the  practical  heat  and  energy  transfer 
systems,  initially  free  from  noncondensable  gases,  would  eventually  be  contaminated  by  them  during  the  course  of 
their  operation. 

Therefore,  we  here  try  to  study  the  effects  of  the  internal  structure  of  the  condensed  phase  on  the  flow  fields  of  a 
mixture,  particularly  the  flow  fields  of  a  condensable  gas  in  the  presence  of  a  small  amount  of  noncondensable  gas. 
Actually,  we  take  up  simple  flow  fields  of  a  mixture  due  to  evaporation  and  condensation  processes  from  or  onto 
the  plane  condensed  phase  with  a  temperature  field  as  its  internal  structure.  For  this  analysis,  the  equation  of  heat 
conduction  with  constant  substance  properties  is  used  for  the  temperature  field  within  the  condensed  phase,  and  for 
the  motions  of  the  mixture  gas  phase  the  Boltzmann  equation  of  BGK  type  proposed  by  Hamel  [4]  is  used  subject  to 
the  diffusive  type  of  boundary  conditions  with  the  parameter  (called  the  imperfectness  parameter  here)  first  introduced 
by  Wortberg  and  his  co-workes  [5]  incorporated  together  with  the  condition  of  the  continuity  of  energy  flow  across  the 
interface. 


GOVERNING  KINETIC  EQUATIONS 


Consider  a  two  phase  system  of  a  gas  phase  and  the  condensed  phase  with  finite  thermal  conductivity.  Let  the  interface 
surface  of  the  condensed  phase  be  at  v  =  0  and  the  other  end  surface  at  x  =  -D.  The  gas  phase,  which  occupies  the 
half- space  (x>0),  consists  of  a  condensable  gas  (A-gas)  and  a  noncondensable  gas  (5-gas).  Initially,  the  condensed 
phase  and  the  mixture  gas  phase  are  in  complete  equilibrium  at  a  temperature  To.  Let  the  pressure  and  the  number 
density  (or  density)  of  the  component  gas  at  this  state  be  Ps0  and  N *  (or  p£),  respectively.  The  superscript  s  denotes 
A-gas  for  s  =  A  and  5-gas  for  s  =  5 .  Suppose  that  at  a  certain  time,  say  t  —  0,  the  temperature  of  the  edge  surface  of 
the  condensed  phase  at  x  =  -D,  which  is  initially  To,  is  suddenly  changed  to  Tc.  The  heat  flow  then  occurs  through 
the  condensed  phase  and  after  a  certain  elapse  of  time  the  temperature  of  the  interface  surface  starts  to  be  changing 
with  time.  This  leads  to  the  onset  of  the  phase  change  process  at  the  interface  of  the  condensed  phase,  giving  then  rise 
to  a  transient  motion  of  the  mixture  gas  phase. 

Now,  for  the  temperature  field  within  the  condensed  phase,  the  equation  of  heat  conduction  with  constant  substance 
properties 


dT 

dt 


d2T 
dx 2 


(1) 


holds  for  —D  <  x  <  0,  where  t  is  the  time;  x  is  the  coordinate;  kc  is  the  thermal  diffusivity  of  the  condensed  phase, 
which  is  here  assumed  to  be  constant  for  simplicity,  and  T  represents  the  temperature  field  within  it.  For  the  description 
of  the  motions  of  the  gas  mixture  in  the  present  problem,  on  the  other  hand,  the  Boltzmann  equation  of  BGK  type  [4] 
is  used,  which  may  be  written  as 
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Ps  = 


ms 

Ms  +  Mt 


(6) 


where  s  appearing  as  superscript  and  subscript  refers  to  A  or  5  and  t  refers  to  5  or  A.  (£*,  fy,  £•)  is  the  molecular 
velocity  vector;  fs,  Ns ,  us ,  Ts  and  Ps  (=  NskTs )  are,  respectively,  the  distribution  function,  the  number  density, 
the  mean  flow  velocity,  the  temperature  and  the  partial  pressure  of  s-component  gas;  Fse  and  FSJ  are  the  local 
Maxwellian  distributions  for  molecules  of  s-component  gas  characterized  by  the  macroscopic  quantities  (Ns ,  us ,  Ts ) 
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and  ( Ns ,  us\  Tst ),  respectively;  ms  is  the  molecular  mass  of  s-gas;  Rs  =  k/ms  is  the  gas  constant  per  unit  mass,  k  being 
the  Boltzmann  constant.  Nskss  and  A/V^  represent  the  number  of  collisions  per  unit  time  of  a  molecule  of  s-gas  due 
to  the  self  and  cross  collisions,  respectively,  thus  Nskss  +  A/V^  being  the  average  collision  frequency  for  a  molecule 
of  s-gas  irrespective  of  its  collision  partners.  The  collision  frequencies  are  related  to  the  transport  coefficients  of  the 
mixture  and  its  component  gases  and,  hence,  kss  and  Kst  (=  Kts),  which  are  assumed  to  be  constant  and  of  the  same 
order  of  magnitude,  can  be  evaluated  from  these  coefficients  at  some  equilibrium  reference  state  (see  e.g.,  Hamel  [4]), 
say  at  infinity,  as 
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where  rjM  and  Dab  are  the  viscosity  and  the  diffusion  coefficients  of  the  mixture,  respectively,  r]s  and  As  being  the 
viscosity  and  thermal  conductivity  of  s-gas. 

Once  the  distribution  function  fs  is  obtained,  the  normal  stress  tsxx ,  the  energy  flux  Es  of  s-gas  can  be  calculated 
from 


Kx  =  -ms  JJJ  ( i' v  -  11s)2  fdt;xdZydt;z 
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(8) 


Now,  as  the  initial  conditions  for  the  present  problem,  we  have 


T  =  T0 


f  = 


Ns 

iVo 


(2nRsT0)3D 


exp 


Zl+Zy+tV 

2  RST0 


for  -  D  <  x  <  0 
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for  the  temperature  of  the  condensed  phase  and  for  the  distribution  function  of  the  component  s-gas,  respectively.  As 
to  the  boundary  conditions,  we  have  on  the  edge  of  the  condensed  phase 


T  =  TC  at  v  =  -D. 


(10) 


At  infinity  (x  — >  oo),  the  condition  for  the  distributionn  function  fs  of  the  component  gas  remains  to  be  the  same  as 
that  given  by  Eq.  (10)  but  for  gx  <  0.  However,  at  the  interface  surface  of  the  condensed  phase  ( v  =  0),  we  have 
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where  N A  and  N ^  represent  the  number  densities  of  the  molecules  of  A-gas  and  B-gas,  respectively,  leaving  the 
interface  surface  after  the  interaction  with  it.  TV  is  the  temperature  of  the  interface  surface  of  the  condensed  phase 
and  is  a  function  of  time  to  be  determined  as  a  part  of  the  solution.  Ac  is  the  thermal  conductivity  of  the  condensed 
phase,  hi  is  the  latent  heat  of  vaporization  per  unit  mass  of  the  condensable  gas  and  c A  is  its  specific  heat  at  constant 
pressure.  The  parameter  ac  represents  a  fraction  of  the  molecules  of  condensable  A-gas  emitted  from  the  interface 
with  an  equilibrium  number  density  whereas  1  —ac  represents  the  remaining  fraction  of  the  molecules  reflected 
from  the  interface  as  if  they  were  noncondensable  gas  molecules.  It  is  a  parameter  associated  with  a  certain  kind 
of  the  imperfectness  of  the  interface  first  introduced  by  Wortberg  and  his  co-workes  [5]  in  their  experimental  and 
approximate  analytical  studies  on  this  kind  of  problems,  and  here  it  is  called  the  imperfectness  parameter.  is  an 
equilibrium  number  density,  which  is  a  unique  function  of  TV,  and  is  given  here  as  the  saturated  vapor  number  density 
corresponding  to  TV  given  by  the  Clapeyron-Clausius  relation  (see  e.g.,  Landau  and  Lifshitz  [6]).  P ^  is  the  saturated 
vapor  pressure  corresponding  to  TV  just  as  P ^  is  to  To.  T  is  a  non-dimensional  parameter,  which  is  here  called  the 
latent  heat  parameter.  The  other  quantity  N^,  which  is  the  number  density  for  molecules  of  A-gas  reflected  from  the 
interface  as  if  they  were  noncondensable  gas  molecules,  is  determined  by  the  condition  of  no  net  mass  flow  across  the 
interface.  N ^  is  also  determined  by  the  condition  of  no  net  mass  flow  of  the  molecules  of  the  noncondensable  5-gas 
across  the  interface.  N ^  and  N ^  are  both  unknown  quantities  to  be  obtained  as  part  of  the  solution. 

Now,  let  us  introduce  here  the  length  scale  L  and  time  scale  for  the  present  problem  defined,  respectively,  by 


L  = 


2  0 


T  /A 

and  V1  =  L  =  0 

0  {2RaTq)1!2  (SRaTo/t: )1/2 


with 


[A  =  ( *RaTq/i t)1/2 
0  N^km+N^ab 


(12) 


where  l ^  is  the  mean  free  path  of  the  molercules  of  the  condensable  A-gas  at  the  initial  state  and,  therefore, 
represents  their  mean  collision  time.  With  these  scales,  the  appropriate  nondimensionalization  of  the  set  of  the 
governing  equations,  initial  and  boundary  conditions  gives  us  the  following  parameters  which  characterize  the  motions 
of  the  mixture  and  its  component  gases: 
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where  TA  and  /cA  are  the  thermal  conductivity  and  the  thermal  diffusivity  of  the  condensable  gas  at  the  initial  state, 
respectively,  k a  being  defined  by  k a  =  AA/(pAcA).  The  specification  of  the  above  parameters  gives  the  flow  field  as 
functions  of  t/rA  and  x/L. 

In  order  to  make  the  numerical  analysis  of  the  present  problem  easier,  the  following  new  functions 


d%yd%z 


and 


"  hs 

Hs 

Hst 


d£yd£z 


(14) 


have  also  been  introduced  here  (see  e.g.,  Chu  [7]).  Rewriting  the  original  equations  with  these  new  variables,  we 
then  obtain  the  simultaneous  equations  for  the  four  unknown  functions  (gA,  hA,  gB,  hB)  of  three  independent  variables 
(t,  v,  gx)  instead  of  having  two  unknown  functions  (/A,  fB)  of  five  independent  variables  (t,  x,  £x,£y,  &)• 


RESULTS  AND  DISCUSSION 


Actual  calculations  have  been  carried  out  based  on  a  simple  difference  scheme  for  more  than  20  cases,  setting 
D/L  =  20.0  and  kaa  =  kab  =  kbb ,  for  simplicity.  Several  of  these  cases  will  be  discussed  here.  For  this  discussion, 
the  fluid  dynamic  quantities  of  the  mixture  5,  A,  p,  E ,  u  and  T,  which  are  defined  by 


P  =  PA+PB,  N  =  Na+Nb,  p=pA+pB,  E  =  Ea  +  Eb,  u=—ua  +  —  ub,  T  =—Ta  +  —  Tb 

y  y  y  p  p  N  N 


(15) 


are  also  introduced.  The  corresponding  quantities  with  suffix  0  represent  the  values  at  the  initial  state  (hereafter,  the 
same  applies  to  the  other  quantities).  cs0  is  the  sound  speed  of  s-gas  at  the  initial  state  defined  by  cs0  =  ( yRsTo )1/2. 
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Incidentally,  the  corresponding  sound  speed  of  the  mixture  Co  is  Co  =  (y^7o)1/2  =  c^(R/Ra)1^2 ,  R  being  the  gas 
constant  per  unit  mass  of  the  mixture  given  by  R  =  (pA  /po)Ra  +  (Pq/Po)Rb • 

One  of  the  typical  transient  flow  fields  due  to  evaporation  process  at  the  interface  surface  is  shown  in  Fig.  1  in 
terms  of  the  temperature  distributions  and  in  Fig.  2  in  terms  of  the  pressures  of  the  mixture  and  its  component  gases. 
Figure  1  shows  how  the  jump  in  temperature  occurring  at  the  interface  varies  with  time  and  how  the  temperature 
distribution  is  being  formed  in  the  flow  field  owing  to  the  propagation  of  the  produced  shock  wave  followed  by  the 
contact  region.  Specifically,  the  right  figure  shows  the  transition  of  the  formation  of  a  shock  wave  with  time  and  at  about 
t/rfi  =  2000  the  shock  wave  almost  fully  establishes  itself,  although  the  rear  portion  of  it,  which  is  still  curved,  is  being 
augmented  gradually  by  the  continuous  overtaking  of  the  compression  waves  produced  at  the  interface  (cf.  Onishi  and 
Yamada  [2]).  The  magnitude  of  the  temperature  jumps  first  increases  from  zero  and  then  approaches  a  constant  value 


FIGURE  1.  Transient  distributions  of  the  temperature  of  the  mixture  due  to  the  evaporation  process  at  the  interface  at  x  =  0. 
TJT0  =  2.0,  D/L  =  20.0,  Ac/Aa  =  50.0,  kc/k$  =  0.5,  T  =  5.0,  NB/NA  =  0.01,  mB/mA  =  2,  kaa=kab  =  kBB  and  ac  =  TO.  The 
numbers  in  the  figures  indicate  the  time  t/r^.  Note  that  the  scales  of  abscissa  of  the  three  figures  are  different.  T /7b|x=o  =  1.0349 
at  t/r*  =  200;  T/T0\x=0  =  1.1146  at  t/rA  =  500;  T/T0\x=0  =  1.1519  at  t/r*  =  1000;  T/T0\x=0  =  1.1584  at  t/r^  =  3000. 


FIGURE  2.  Transient  distributions  of  the  pressures  of  the  mixture  and  its  component  gases  due  to  the  evaporation  process  at  the 
interface  at  x  =  0.  Tc/T0  =  2.0,  D/L  =  20.0,  Ac/Aa  =  50.0,  kc/k£  =  0.5,  T  =  5.0,  NB/NA  =  0.01,  mB/mA  =  2,  km  =  km  =  kbb 
and  ac  =  1.0.  The  numbers  in  the  figures  indicate  the  time  t/r^.  Pa/Pq\x=o  =  2.0599,  Pb/Pq\x=o  -  0.0,  P/Po\x=o  =  2.0395  at 
1/^  =  3000. 

with  time,  which  means  that  the  flow  field  has  approached  a  final  state  (not  a  steady  state,  strictly  speaking).  In  the 
figures,  the  steep  gradients  of  the  temperatures  and  pressures  can  be  recognized  in  a  small  region  (x/L  ~  0(1))  in 
the  immediate  vicinity  of  the  interface.  This  region  is  the  so-called  Knudsen  layer,  the  thickness  of  which  is  of  the 
order  of  the  molecular  mean  free  path,  and  it  always  manifests  itself  in  the  close  vicinity  of  the  boundary  surface.  In 
the  uniform  region  (0  <  x/L  ~  600  at  t/rA  =  3000),  which  connects  the  contact  region  and  the  Knudsen  layer,  the 
noncondensable  gas  contained  in  a  small  amount  initially  is  almost  zero  because  it  is  driven  almost  completely  into 
the  contact  region  owing  to  the  evaporation  flow.  The  roles  of  the  latent  heat  parameter  T  and  the  thermal  conductivity 
ratio  Ad Aa  are  qualitatively  the  same  in  the  present  mixture  case  as  in  the  pure-vapor  case  studied  earlier  by  Onishi 
and  Yamada  [2],  i.e.,  the  larger  the  values  of  T,  the  weaker  become  the  phase  change  processes  and  hence  smaller 
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become  the  mass  and  energy  flows  and  the  larger  the  thermal  conductivity  ratio  Ac/A^ ,  the  stronger  become  the  phase 
change  processes  (the  physical  explanation  is  given  there).  However,  when  T  becomes  smaller  than  a  certain  definite 


FIGURE  3.  The  total  mass  flow  pu/poc^  and  energy  flow  E*  =  E / Pq{2RaTq)1!2  of  the  mixture  at  the  steady  state  versus  the 
latent  heat  parameter  T.  Tc/Tq  =  2.0,  D/L  =  20.0,  N^/Nq  =  0.01,  m^/m^  =  2.0,  kaa  =  kab  =  kbb  and  ac  =  1.0.  The  white  circle: 
Ac/Aq  =  50.0,  kc/k£  =  0.50;  The  balck  circle:  Ac  /Aq  =  30.0,  kc/k^  =  0.32.  The  existence  of  the  maxima  of  the  mass  and  energy 
flows  can  be  expected.  The  graph  on  the  flush  right  shows  the  temporary  negative  mass  flow ,  the  parameters  for  which  are  the  same 
as  in  Figs.  1  and  2  except  that  here  T  =  0.01.  The  values  of  this  negative  mass  flow  are  pu/poc^\x=o  =  -0.0027  at  t/r^  =  300  and 
pu/poc^\x=o  =  -0.0021  at  t/r^  =  500. 

value,  the  mass  and  energy  flows  decrease  again,  which  means  that  there  exist  their  maximum  values  at  a  certain  value 
of  T.  This  behavior  is  shown  in  Fig.  3.  This  is  a  striking  contrast  to  what  we  have  studied  so  far  in  which  the  condensed 
pahses  have  no  internal  structures  at  all. 

Finally,  the  negative  mass  flow ,  another  feature  of  the  present  problem,  manifests  itself  during  a  certain  period  of 
time  when  the  latent  heat  parameter  T  becomes  extremely  small  but  it  disappears  soon  as  time  goes  on.  A  typical 
example  of  this  is  shown  in  the  graph  on  the  flush  right  of  Fig.  3.  The  existence  of  this  negative  mass  flow  at  the  steady 
state  has  first  been  discussed  analytically  by  Sone  &  Onishi  [8].  The  explanation  is  also  given  in  Onishi  &  Ooshida  [9]. 
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